Graphene Antidot Lattices — Designed Defects and Spin Qubits 
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Antidot lattices, denned on a two-dimensional electron gas at a semiconductor heterostructure, are 
a well- studied class of man-made structures with intriguing physical properties. We point out that 
a closely related system, graphene sheets with regularly spaced holes ("antidots"), should display 
similar phenomenology, but within a much more favorable energy scale, a consequence of the Dirac 
fermion nature of the states around the Fermi level. Further, by leaving out some of the holes one 
can create defect states, or pairs of coupled defect states, which can function as hosts for electron 
spin qubits. We present a detailed study of the energetics of periodic graphene antidot lattices, 
analyze the level structure of a single defect, calculate the exchange coupling between a pair of spin 
qubits, and identify possible avenues for further developments. 
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Graphene is the rapidly rising star of low-dimensional 
materials. Following the initial reports on fabrication 
by mechanical peeling [1] and epitaxial growth [2], this 
exceptional material has stimulated considerable experi- 
mental and theoretical research [4] as well as proposals 
for novel electronic devices [5]. The promising prospects 
for graphene devices are based on several remarkable 
properties. Mainly, the sample quality and mobility (ex- 
ceeding 15000 cm 2 /Vs [3]) can be very high. In addition, 
patterning of such monolayer films by e-beam lithogra- 
phy [E B] with features as small as 10 nm 0, Q is pos- 
sible. Very recently, spintronics devices have been con- 
sidered 0. The incentive for graphene based spintronics 
lies partly in the long spin coherence time that is charac- 
teristic of carbon-based materials. This also has obvious 
advantages within the field of solid-state quantum infor- 
mation processing, where confined electron spins have 
been promoted as carriers of quantum information [9[. 
Being a light element, carbon has a rather small spin- 
orbit coupling, and, moreover, the predominant 12 C iso- 
tope has a vanishing hyperfine interaction. This makes 
graphene, at least in principle, a superior material com- 
pared to existing quantum computing implementations 
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Antidot lattices, defined on semiconductor het- 
erostructures, display many intricate transport proper- 
ties, in particular in magnetic fields where the compet- 
ing length scales lead to rich physics [l2[. In this Letter 
we wish to draw attention to the possibility of forming 
antidot lattices on graphene. As mentioned above, state- 
of-the-art e-beam lithography has been used to carve 
graphene nanoribbons with feature sizes down to tens of 
nanometers. We propose to use similar techniques to cre- 
ate regular holes in the graphene sheet, in order to form 
antidot lattices. The antidot lattice has the important 




JUU% w wrk 
.+.*.*.« i , . *■*»%*. 

'*:*;# 1 055* 

•******i»i*.*«* 



•*•*;* *;*;*; 
****** i 

**•*•• »***•* 

xosccr 



msssr 



FIG. 1: Illustration of the triangular antidot lattice (a) with 
a unit cell characterized by side length L and hole radius 
R (b). In (c), several examples with corresponding {L, R} 
parameters are shown. 



consequence that it turns the semi-metallic graphene into 
a gapped semiconductor, where the size of the gap can be 
tuned via the antidot lattice parameters. As our analysis 
shall show, this electronic structure can be manipulated 
further so as to create coupled electron spin qubits, thus 
suggesting that these perforated graphene sheets are a 
promising platform for a large-scale spin qubit architec- 
ture. Localized spin qubit states can be formed in the 
antidot lattice by deliberately omitting some of the an- 
tidots. This idea has previously been analyzed for the 
two-dimensional electron gas in e.g. GaAs heterostruc- 
tures [13[. As we will now argue, moving to graphene 
has three major advantages: (i) increased coherence time; 
(ii) favorable energy scale of the defect states; and (iii) 
increased lateral confinement. 

The proposed antidot lattice is simply a triangular ar- 
ray of holes in a graphene sheet, as illustrated in Fig. 
[T^l. The lattice consists of hexagonal unit cells as shown 
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in Fig. [TJd, in which a roughly circular hole is created. 
We characterize the structure by the side length L of 
the hexagonal unit cell and the radius R of the hole, 
both measured in units of the graphene lattice constant 
a « 2.46 A. A lattice is designated by the notation 
{L,R}. Note that while L is an integer, R can be non- 
integer. As is evident from the examples in Fig. [Tfc, L is 
equal to the number of carbon atoms in the outermost 
row of the hexagon. Also of importance are the total 
number of sites in the unit cell A^otai (equal to the num- 
ber of atoms before the hole is made) and the number of 
removed atoms iVR e moved- As an example, for the {7,3} 
lattice AT Total = 294 and A/R emove d = 60. Below, results 
for structures with L < 14 and varying R have been com- 
piled taking care that no dangling bonds are formed, i.e. 
that all atoms have at least two neighbors. While these 
structures are too small for present-day lithography, re- 
sults for realistic structures are easily obtained by simple 
scaling laws, as demonstrated below. 

We model the structures using a tight-binding (TB) de- 
scription considering a single 7r-orbital on each site and 
assuming a nearest-neighbor hopping integral of — /?, with 
(3 « 3.033 eV [14j. In this description, energy levels are 
always distributed symmetrically above and below zero, 
which defines the Fermi level in the undoped case. The 
TB approximation is necessary due to the large antidot 
cells. It is known to accurately reproduce the low-energy 
part of the density-functional (DFT) band structure of 
graphene [15[. Edges, however, require a modification of 
hopping integrals near the edge to ensure agreement be- 
tween DFT and TB calculations [la]. We have checked 
that the computed band structures are generally robust 
against such modifications, which simply produce a mi- 
nor additional opening of the band gap. The electronic 
band structure and density of states for the {7, 3} struc- 
ture are illustrated in Fig. O Importantly, a substantial 
energy gap of approximately 0.73 eV opens around the 
Fermi level [ljj. Hence, as hinted above, the periodic per- 
turbation turns the semi-metal into a semiconductor. In 
the top panel of Fig.[3j band gaps E g of several structures 
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are plotted versus the quantity Af R ^ moved /iVrotai- When 
plotted in this manner, a roughly linear behavior is ob- 
served. This simple result may be rationalized within 
the linearized Hamiltonian approximation treating elec- 
trons as massless Dirac fermions subject to the periodic 
perturbation of the antidot lattice. In this description, 
the wave function is a two-component spinor represent- 
ing the two sublattices. The corresponding Hamiltonian 
is the 2x2 matrix operator 

H = ( V ( X > V) V F(Px - iPy)\ 

\v F (px + ip y ) V(x,y) J 

where V is the periodic antidot potential, p is the 
momentum operator, and the Fermi velocity vp = 
\ /r 3(3a/(2H) « 10 6 m/s. In the absence of a potential, 
the energy eigenvalues are simply E = ±hvp\k\. If the 




FIG. 2: Energy band structure and associated density of 
states for a {7, 3} antidot lattice. The notation T, M, and 
K refers to high symmetry points of the Brillouin zone. 

potential is approximated by infinite barriers at the po- 
sitions of the antidots, the eigenvalue problem is reduced 
to the form 

v 2 f (pI+P 2 v W = E 2 iP, (2) 

with the boundary condition that ip vanishes in the bar- 
rier region. The equation is mathematically similar to 
the usual effective mass equation. For an antidot lat- 
tice in a usual semiconductor material such as GaAs, 
simple scaling arguments lead to a band gap varying 
as Eg oc ^^^/(ARemovedMTotai), where Axotai is the 
area of the unit cell and AR emove d is the area removed 
inside each unit cell. In graphene, a similar behav- 
ior is expected except that the linear band structure 
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changes the prefactor from ^ Xotal to ^4 Total , i.e. E g cx 

^Total# (^Removed / ^Total ) OC ^Total # (^Removed / ^Total ) • 

The fit in Fig. [3] shows that g approximately fol- 
lo ws a square ro ot behavior g (N Removed /AT Total ) ex 
V^Removed/^Totai- Thus, the net result is a gap varying 

1/2 

as Eg w K x ^Removed/A^Totai with a constant K w 25 

1/2 

eV. For large unit cells, N^ emoved A^Totai is small and in 
this case the linear fit is an excellent approximation. The 
weaker scaling (^xotai i ns tead of A^otai) of graphene is 
very favorable for the purpose of obtaining large band 
gaps even for relatively large structures. The practical 
limits of present day e-beam lithography probably re- 
strict the obtainable size of the unit cell to around 10 nm 
across corresponding to a total number of carbon atoms 

Of A^Total ~ 3000. Assuming ^Removed ~ N To tal/4 We 

find a substantial gap of 0.23 eV. Hence, band gaps much 
larger than the thermal energy at room temperature are 
certainly realistic. This feature, which is a direct con- 
sequence of the massless Dirac fermion behavior, is very 
important for the feasibility of the graphene based de- 
vices considered here. 

We now turn to the role of intentional defects in the an- 
tidot lattice produced by leaving one or several unit cells 
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FIG. 3: Compilation of energy gaps (upper panel) and defect 
state binding energies (lower panel). When displayed versus 
^Removed /-^Totai? very simple scaling is observed. Note that 
^R^Loved/^Totai is small for realistic structures. 



intact, i.e. without a hole. Such defects may support lo- 
calized electronic states and may consequently be utilized 
for electron spin qubits, as we will now demonstrate. An 
example of single and double defects for the {5,2} struc- 
ture is shown in Fig. [U For isolated single defects, we 
compute localized states by periodically replicating the 
super cell consisting of one intact and six perforated cells 
illustrated in the figure. The states are sufficiently lo- 
calized that cross-talk between neighboring super cells 
is negligible. Periodicity is not crucial for the appear- 
ance of bound states [13[ . Defect states are identified by 
an energy lying in the fundamental energy gap, i.e. the 
gap containing the Fermi energy. In fact, other energy 
gaps may exist as illustrated in Fig. [2j here we focus 
solely on states in the fundamental gap. If the gap is 
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sufficiently large (i.e. if ^V Removed /^Totai is large) several 
defect states are supported. In the lower panel of Fig. O 
a compilation of binding energies for the three lowest de- 
fect states is shown. We define the binding energy £?Bind 
as the downwards shift of the defect state energy mea- 
sured from the conduction band edge. Hence, a defect 
state at the Fermi energy would have a binding energy 
of Eg/2. For small band gaps, only a single defect state 
is supported but several defect states appear in an irreg- 
ular pattern as the confinement increases. Note that the 
scatter in the data points in the plot reflects actual vari- 
ations and not computational inaccuracy. Importantly, 
the binding energy in the limit of small band gaps is seen 
to approach a constant fraction ~ 0.07 of the energy gap. 
Hence, for the 10 nm unit cell considered above, a defect 



FIG. 4: Single (left) and double (right) defects for the {5, 2} 
antidot lattice. To compute defect states, super cells con- 
taining defects surrounded by six intact units are repeated 
periodically. 



state would be bound by roughly 16 meV. This implies 
that liquid nitrogen cooling should be sufficient to ob- 
serve these states. 

Next, we consider two tunnel coupled defect states 
in a "double defect", illustrated in Fig. [H With an 
electron occupying a non-degenerate state in each de- 
fect, the spins of the two electrons couple due to the 
exchange interaction JSi • S2. If the two single-defect 
states are energetically aligned, the exchange coupling is 
given as J = 4t 2 /U according to the Hubbard approxi- 
mation. Here, t is the tunnel coupling between the two 
defect states, and U is the single-defect Coulomb energy. 
As discussed in Ref. [9], the exchange coupling consti- 
tutes a key element in quantum computing architectures 
based on electron spins as qubits, enabling interactions 
between different qubits. Importantly, the exchange cou- 
pling can be controlled with external gate potentials. 
Metallic gates could be realized by lithographic methods 
and placed either below or on top of the graphene sheet 
but will not be considered further here. For evaluation 
of the exchange coupling, we calculate the single-defect 
Coulomb energy U by the method presented in Ref. 18 
(ignoring overlap between different atomic 7r-orbitals) us- 
ing the Ohno form to interpolate between the intra- and 
long-range inter-atomic Coulomb coupling. A Hubbard 
U n for carbon 7r-orbitals of 20.08 eV [18[ and dielectric 
constant of 2.5 [19[ (as appropriate for graphene on Si02) 
are applied. The tunnel coupling t is extracted from the 
single-particle energy spectrum. 

Our findings for the Coulomb energy U are illustrated 
in Fig. [5l In the plot, Rd is the effective defect ra- 
dius calculated by including half the area of the sur- 
rounding cells and writing the total area as ttR 2 d . The 
smallest C/'s are found for the least localized states for 
which U scales as the expected R^ 1 . The inset shows, as 
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FIG. 5: Compilation of Coulomb energies U for single defect 
states showing roughly linear scaling with R^ 1 for delocalized 
states. Inset: Level structure for single- and double defects in 
a {12,7} lattice. The Coulomb energy for this case is indicated 
by the circle. 



an example, the single-electron level diagram for single- 
and double defects in a {12,7} lattice. This structure 
has 7Vxotai= 864 and A^R emove d = 348 and supports two 
single-defect states. Of these, the upper one is non- 
degenerate and the Coulomb energy is 0.315 eV. Due to 
the large double defect super cell, this is about the largest 
structure that we have been able to analyze. As shown, 
the level splitting corresponds to a tunnel coupling of 
t ~ 2 meV between the two non-degenerate single-defect 
states. Hence, based on the {12,7} values we may esti- 
mate the exchange coupling to be on the order of J ~ 50 
fieV. Naturally, this value could be tuned by appropri- 
ate design of the barrier region that, for simplicity, has 
been constructed from two intact unit cells. Also, going 
to larger single defects would decrease U and, in turn, 
increase J. Note, however, that t depends exponentially 
on barrier width whereas U is only weakly dependent on 
geometry. Hence, the geometric influence on J will be 
determined mainly through t rather than U. 

We believe that the approach outlined above can be 
extended to more complicated structures. Going from a 
single pair of spin qubits in an isolated double defect to 
several coupled spins could be achieved with little added 
complication. Similarly, a double defect could be re- 
placed by a linear array of defects. Hence, the number of 
qubits can be increased essentially without complicating 
the fabrication procedure. In practice, excellent control 
of the e-beam lithography process remains a critical issue. 

In summary, we have shown that antidot lattices pave 
the way for controlled manipulation of the electronic 
properties of graphene sheets. The material can be ren- 
dered semiconducting with a significant and controllable 
energy gap. The magnitude of the gap is explained by a 



simple scaling argument and could reach several tenths of 
eVs for realistic structures. Introducing defects into the 
antidot lattice leads to the formation of localized elec- 
tronic states. Combined with the extremely long spin 
coherence time of carbon-based materials this could lead 
to a practical realization of spin qubits. With a properly 
designed double defect, two-electron states derived from 
defect levels near the Fermi level are found to fulfil the 
requirements for such qubits. 
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